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ABSTRACT 






The purpose of this thesis is to find an automated way to identify push-to-talk 
transmitters using a wavelet or a spectral correlation based approach. In the Wavelet 
approach, a distance algorithm is applied to the wavelet scales of the signal and the 
template. One signal from each transmitter signal set is taken as a template. The 
distance algorithm computes the distance between the local extrema of the wavelet 
coefficients of the template and the signal. Results show that the Wavelet Transform 
(WT) distance algorithm is able to classify the four signal sets accurately. Good 
identification results are achieved even at low signal-to-noise ratios. In the spectral 
correlation approach an averaged template for each signal set is used. The spectral co- 
efficients for templates and signals are computed by extracting the magnitude squared 
of the Fast Fourier Transform (FFT) of the data. This method performs better for 
most signals than the wavelet method because it can identify at lower signal to noise 
levels than the wavelet method does. 
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I. INTRODUCTION 



Some characteristics of signals generated by a transmitter are believed to be 
unique to a given transmitter. The transient response is defined as the change in car- 
rier strength from the off-state to the on-state, and vice versa. The turn-on response 
(from the off- state to the on-state) is unique for every transmitter, regardless of its 
make or model. By analyzing the turn-on response we can map a given intercepted 
signal to a particular transmitter. Therefore, if we can reliably extract the unique 
parameters of the turn-on transient, then we can identify a given transmitter. 

Time-frequency analysis of stationary signals is a well- researched subject. This 
type of analysis uses the Fourier Transform (FT) method which is appropriate for 
stationary signals. However, since the FT method uses a complex exponential ba- 
sis function that exists over infinite time, we cannot use the FT method for non- 
stationary signals. In 1946, Gabor [1] introduced the sliding time window (Gaussian) 
to gain time information from the FT method. This modified FT method is called 
the Short-Time Fourier Transform (STFT). It uses a modulated complex exponential 
as a basis function and processes only the windowed data segment. This method, 
however, requires the signal to be stationary within these small portions, otherwise it 
will have limitations in reflecting the time-evolution of frequencies within the window. 

The Wavelet Transform (WT) method is a more suitable technique for transient 
signal analysis than the STFT method in the sense that it provides more conclusive 
information regarding time evolution. Basis functions of the WT, unlike the FT, have 
a very short duration and are nonzero only for that short duration. This compact 
support makes the WT localized in frequency and time. Moreover, if we use wavelets 
we can choose the particular wavelet function from a large number of compactly 
supported wavelets that can be used as orthogonal basis functions. 
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Spectral correlation analysis provides a useful tool for the classification of tran- 
sient signals (i.e., finite energy). Spectral coefficients are obtained by using the Fast 
Fourier Transform (FFT) of the data and extracting the magnitude square of the 
coefficients at the location of interest, as well as the spectral location. This technique 
correlates spectral coefficients of signals with those obtained from selected templates. 

This thesis consists of five chapters. Chapter I introduces the topic. Chapter II 
explains briefly wavelet analysis, the processing of the wavelet transform output, and 
the spectral cross-correlation technique. Chapter III explains the signal precondition- 
ing. Chapter IV details the processing results and Chapter V presents conclusions 
and recommendations for future studies. 
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II. TRANSIENT SIGNAL ANALYSIS 



A. WAVELET TRANSFORM (WT) 

The Wavelet Transform (WT) is an alternative technique to the classical Short 
Time Fourier Transform (STFT) (i.e., the Gabor Transform [1]) which is used for 
the time-frequency analysis of non-stationary signals [2]. The STFT uses a single 
fixed analysis window, while the wavelet transform uses various analysis windows in 
the time domain. The wavelet transform bandwidth is inversely proportional to the 
width of the time window. Hence, the more narrow the wavelet function window 
width, the higher the center frequency is and the wider the bandwidth is. This is the 
reason why the WT is known as “constant Q” filtering [3]. The wavelet transform 
uses dilations and translations of a basis function known as the mother wavelet. The 
family derived from this basis function is generally of the form: 

Mt) = ^ (^) - (i) 

where if) is the mother wavelet, a is the scaling parameter, and b is the shifting or 
translation parameter. We numbered the octaves such that as the frequency increases 
the octave increases too. The constant l/y/a is used for energy normalization. The 
family of the basis functions include short duration/high frequency functions and 
long duration/low frequency functions. The dilations and compression of the wavelet 
function ?/>(£) enable the tiling of the time-frequency plane, as shown in Figure 1. 

1. The Continuous Wavelet Transform 

The Continuous Wavelet Transform (CWT) can be defined as: 

CWT x (a, b) = 4= r 0, , b (t)x(t)dt , (2) 

yJCL J —oo 

where il> a b(t) is given by Equation (1). 
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frequency 




time 

Figure 1: Tiling of the time-frequency plane. 



A function i/>(t) can be defined as a mother wavelet only if it is oscillatory, 
decays to zero at the end points, and has a zero mean [4]. The inverse of the Wavelet 
Transform is given by 



*W = 7 hr CWT x (a,b)tl> a „(t) 

J—oo 



where C ^ is given by 






dadb 
2 5 



= /° 



n 



dfl , 



( 3 ) 



( 4 ) 



and is the Fourier transform of Equation (3) is valid only if is an 

admissible wavelet function which requires that < oo. This implies that ip a b{t) has 
a zero mean. The function ip a b{t) has compact support, hence it acts as the impulse 
response of a bandpass filter that decays very fast to zero [2]. 

2. The Discrete Wavelet Transform 

From Equation (1), the Discrete Wavelet Transform (DWT) can be obtained 
by sampling the parameters a, b at discrete values which are given as 
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b — kbo&Q^ , 



(5) 



where j is the octave number, k is the sample index, j,k 6 Z and a 0 , 6 0 are pre-selected 
constants. We can express the discrete wavelet functions as 

*l>jk(t) = a 3 J\[a 3 0 t “ kb o) i (6) 

where j , k E Z. The signal can be expressed as a function of this wavelet basis: 

= ( 7 ) 

3 k 

where djk is the wavelet coefficient defined by: 
djk = (^(Oj 

= a J J 2 [ x(t)ip*(a J 0 t — kb 0 ) • (8) 

J — oc 

Under specific conditions detailed in [5], the set of ipjk(t) is considered as an orthonor- 
mal set of basis for L 2 (R). Therefore, any function x(t) in L 2 (R) can be represented 
as in Equation (7). 

The uncertainty principal states that, for a given transform pair, x(t) X(lj ), 

the inequality below holds 

w* > i , (9) 

where a t and <j w are defined as [4]: 

mmt . jw’ixmi , im 

' /ix(oi 2 d< ’ - i\x(u)\^ ■ 

We note the integration equation limits ioo are suppressed. As seen from Equation 
(9), there is a lower bound for the product cr f (7 w , hence if the signal is short in time, 
it will be wide in frequency. This means that we cannot find a wavelet function or 
Fourier basis that has an arbitrary time width and frequency bandwidth, i.e., even 
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for the wavelet transform analysis the time localization is a function of the frequency 
resolution. Accordingly, the high octave corresponds to a narrower wavelet with 
better time resolution than at a lower octave. 

3. Multiresolution Signal Analysis 
a. Multiresolution Spaces 

Multiresolution analysis consists of a sequence of closed subspaces [10] 

• • • V-2 C V~! C Vo C Vi C V 2 • • • (11) 

that have the following properties [6] 

• Completeness: 

n Vj = {0} , u Vj = L*(R) jez , (12) 

where D is the intersection of subspaces and U is the union of subspaces. 

• Scaling Property: 

x{t) 6 Vj x(2t ) £ Vj. f_i . (13) 

• Existence of Basis: A scaling function 4> € Vj exists such that the set 

{<M*) = 2 J/ V(2 j t - k)} , Vjez (14) 

forms an orthonormal basis for Vj. Equation (11) is demonstrated via the Venn 
diagram in Figure 2, as the basic idea of multiresolution analysis is that of 
successive approximations. An equivalent way of representing L 2 (R) is shown 
in Figure 3. 

The space IFj, obtained for octave j, defined to be the orthogonal com- 
plement of the space Vj with respect to Vj+ 1, can be written as 

Vj = Vj_i © Wj-i 

Vj-i±Wj-i , 
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(15) 

(16) 




Figure 2: Equation (11) illustrated via a Venn diagram [7]. 
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Figure 3: Multiresolution representation of L 2 (R) [8]. 
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where the symbol, 0, stands for the direct sum, which means that each element of 
Vj + 1 can be written as the sum of the elements of Vj and Wj. Thus, Wj contains the 
“detail” necessary to go from Vj to Vj+i- 

b. Scaling Functions and Wavelets 

The scaling function 4>jk(f) used in Equation (14) is a lowpass function 
because it is used to go from a higher octave number to a lower octave number. 
It spans the coarse approximation space Vj. The wavelet function *l>jk(t) given in 
Equation (6), which for all values of spans the difference information space Wj , 
can be considered as a highpass function. Scaling and wavelet functions can represent 
any signal in Vj+\. Scaling a wavelet simply means stretching (or compressing) it. 
As per Equations (6) and (14), {ipjk(t)} and {^jfc(O) are se ^ s °f translations of the 
dilated functions ipj(t) and (/>j(t), where the scale factor is denoted by the letter j. 
When j > 0, <t>jk(t) compresses and V 3 stretches accordingly, which results in finer 
details. When j < 0, </>jk{t) stretches and Vj compresses, which results in coarser 
details. 

Let <j>(t ), as used in Equation (14), be the scaling function, then its transla- 
tion <t>(t — k) spans V 0 . The scaling function (/>(t) can be written as a linear combination 
of the translations of <f)(2t) as 

*(0 = £M*W2 <-*), kez, (17) 

k 

where h 0 (k) are the scaling function coefficients. 

Let ^(0? as defined by Equation (6), be an element of the subspace Wo. 
Wq itself is a subspace of Vi spanned by <j>(2t). Hence, ^(t) can be written as the 
linear combination of the translation of <j>(2 1) as 

r/>(t) = J2 hi(k)<i>(2t — k ) , k 6 Z , 

k 
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( 18 ) 



where the h\(k) are the wavelet coefficients. 

Since (f)(t ) and ip(t) span two orthogonal spaces, h 0 (k) and hi(k) should 
be orthogonal. Hence, orthogonality of V 0 and W 0 can be accomplished by requiring 
that 

(h 0 (k), hi(k)) = 0 . (19) 

Any x(t) in L 2 (R) can be represented as 



x (t)= J2 a ^k{i) + Y, J2 

k= — oo j = l k= — oo 



( 20 ) 



with ak and djk being discrete wavelet transform coefficients which are defined by 

/ oo 

x(t)(j> k {t)dt 

-oo 

/ oo 

x(t)rl>jk(t)dt , (21) 

-OO 

where and *l>jk(t) are real functions. According to Equation (20), any function in 
L 2 (R) can be written as a linear combination of the scaling function at a fixed scale 
plus a linear combination of wavelet functions at a higher scale. For more details 
regarding the properties of wavelet transform, see Payal [8]. 

A typical wavelet set up is shown in Figure 4. This is according to the 
following equations [9]: 



CLjk — ) ^ A-o ( 77T. 2 k ) dj _|_ i j 

m 

djk = "y ", hi(m 2k)aj + i ir 

771 

where h 0 (k) and hi(k) are given by 

h 0 (k) = h 0 (-k) 

hi(k) = hi(—k ) . 



(22) 

(23) 



(24) 
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where h 0 (k) and hi(k) are the impulse response time reversals. The wavelet imple- 
mentation for three levels is shown in Figure 5. The notation LP stands for lowpass 
FIR filter whose weights are /i 0 (— &), and HP stands for highpass FIR filter whose 
weights are hi(—k). 




Figure 4: One stage of multiresolution signal decomposition. 




Figure 5: Three level multiresolution analysis. 
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4. Processing of Wavelet Transform Output 

a. Reduced Set Representation 

The method introduced previously in the preprocessing of the wavelet 
transform output [8] is applied in this thesis, where the extrema of the wavelet coeffi- 
cient are used to identify the original signal. The wavelet coefficients are replaced by 
their extrema and the Euclidian distance between these extrema is measured. The 
computer programs related to preprocessing, wavelet processing, and cross-correlation 
of the spectral coefficient are presented in Appendix A. 

b. Ranking and Matching Algorithm 

Once the peaks, which are the extremal values of the wavelet coefficients, 
are provided for the templates and the signals, they are ranked by their amplitudes. 
Then, matching the ranked peaks to form pairs, the peak with the highest rank in 
one set is matched to the peak with the highest rank in the other. Finally, the next 
in rank is matched by the next of the other set and so on. 

c. Distance Measure 

The distance measure, which was introduced in 1992 by Aware, Inc., is 
used in a modified form for the classification of signals in this thesis work. The 
Euclidian distance is defined as 



d(a’,V)= Y, [ W b (4-^)1 ^ (25) 

(k,£) are locations of the matched peaks 



where a 3 and b 3 are wavelet maxima at scale j. W kt is the weighting factor at scale 
j for the relative distance between the similar coordinates of the matched peaks, and 
a k and b J k are the values a 3 and b 3 at momentary location k and t. The weighting 
factor W k l is defined as 



wi, t = \ |ni 7”*i ; 



k±l 

k = £ ’ 



(26) 
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B. SPECTRAL CROSS-CORRELATION 



1. Introduction 

Spectral correlation analysis provides a useful tool for the classification of 
transient signals. These signals (data) are functions of time and are given in digitized 
form. Spectral coefficients are obtained by applying the Fast Fourier Transform (FFT) 
to the data, and extracting the magnitude squared of the coefficients at the location 
of interest, as well as the spectral location [10]. 

Each of the digitized transient signal is 4096 data points long. When the 
spectral coefficients were plotted it was found that all data of interest was in the first 
32 spectral locations. The plots were inspected to confirm that none of the remaining 
locations contained valuable data. To be considered valuable, the spectral density 
must exceed 10% of the peak value. The resultant is plotted for all data sets of all 
transmitters (see Figure 6). We note the relevant information always exists only in 
the low frequency region. 

For this analysis, an averaged template for each signals set is obtained by 
averaging all spectral densities of a given transmitter (see Figure 7). Then the spectral 
coefficients of all signals of interest are obtained. These spectral coefficients are 
examined to determine whether or not this signal belongs to the set which created 
the template. The first step in the processing chain is to find the power in the first 
32 spectral locations given by 



A(k) = 



4095 kn 

E«.W^ 

n=0 iV 



for k = 0, 1, • • • , 31 ; 



(27) 



where a s (n) is the discrete time signal of the data (signal) and N is the data equal to 
length 4096. The spectral coefficients are computed for all templates and signals of 
interest. The degree of similarity between two vectors A te mpiate(k) and A s i gna ](/:) can 
be determined by correlating the vectors under consideration. 
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frequency (K) 

(a) 



x 10 4 




(c) 




frequency (K) 

(b) 




frequency (K) 

(d) 



Figure 6: Superimposed spectra for transients contained in (a) Transmitter 
1, (b) Transmitter 2, (c) Transmitter 3, (d) Transmitter 4. 
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Figure 7: Averaged spectral plot for the four sets of transients: (a) average 
template for TR1, (b) average template for TR2, (c) average template for 
TR3, (d) average template for TR4. 
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2. Normalized Cross-Correlation Coefficient 

Prior to cross- correlating the two vectors A t (k) and A s (k ), the mean value of 
each vector (i.e., template and signal) is removed. The coefficient obtained is the 
normalized correlation coefficient given by 

EMk)A s (k) 

Pxy = 3i t=0 . (28) 

, £K(*)| 2 f> s (i :)| 2 

\ k = 0 k = 0 

Thus, the range of this normalization will keep the cross- correlation coefficient be- 
tween — 1 and +1 and will aid in an automated decision. If p, the normalized cross 
correlation coefficient is one it indicates that the two signals are identical (even if they 
are scaled versions of each other). A high value implies that the template and the 
signal are from the same set, while a small value implies that they belong to different 
sets. 
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III. SIGNAL PRECONDITIONING 



A. TURN-ON AND TURN-OFF SIGNALS OF PUSH-TO-TALK 
RADIOS 

The data was collected and recorded by the Naval Security Group Activity, 
Charleston, SC. Nine samples of turn-on and turn-off transients of four different trans- 
mitters were collected. The recordings from each transmitter are numbered from one 
to nine. Figure 8 shows the first turn-on transient from each transmitter. The signals 
were collected by an antenna, processed with a radio receiver of a carrier frequency 
equal to 138.525 MHz, filtered with a 1 MHz bandwidth filter (BW), and digitized 
with a sampling frequency of 5 MHz at a center frequency of 1.075 MHz. Figure 9 
shows the first turn-off transients from each transmitter. Significant information for 
the identification can be found in the envelope, which contains low frequency com- 
ponents. Wavelet transforms are not practical for the analysis of low frequency, but 
are well- suited for short duration phenomena. Accordingly, the data is transformed 
into a form suitable for wavelet transforms. This is done in the preprocessing phase. 
It should also be noted that the signal-to-noise ratios of the signal recorded are not 
known, and a denoising process would enhance the identification performance. All 
the radios are Motorola radios. Each is identified by its name, number, and model as 
tabulated in Table 1. 

B. PREPROCESSING PHASE 

The signals go through four stages in the preprocessing phase. The four stages are: 
taking the envelope, filtering, differentiating, and a final filtering. Before taking the 
envelope, the D.C. term is removed. The first filter, which filters the envelope, is a 
100-point boxcar averager. The second filter, which is used after differentiation, is a 
50-point boxcar averager. The filters sizes were experimentally determined [10]. The 
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Amplitude Amplitude 



processing of the data deals only with the turn-on transient because the turn-off 
transient does not permit identification. Figure 10 shows the first preprocessed signal 
of each of the four transmitters. All of the four final pulse-shaped signal waveforms 
seem to be suitable for wavelet transform analysis. 




0 2000 4000 6000 0 2000 4000 6000 

K K 

(C) (d) 

Figure 8: Turn-on signal samples from four different transmitters: (a) 

Signal 1 from Transmitter 1, (b) Signal 1 from Transmitter 2, (c) Signal 1 
from Transmitter 3, (d) Signal 1 from Transmitter 4. 
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Amplitude Amplitude 




0 2000 4000 6000 0 1000 2000 3000 4000 

K K 

(c) (d) 



Figure 9: Turn-off signal samples from four different transmitters: (a) 
Signal 1 from Transmitter 1, (b) Signal 1 from Transmitter 2, (c) Signal 1 
from Transmitter 3, (d) Signal 1 from Transmitter 4. 
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Table 1: Name and Model of the Four Transmitters 



Radio 


Model 


Transmitter 1 (Trl) 


Maxtrac 


Transmitter 2 (Tr2) 


Saber 


Transmitter 3 (Tr3) 


HT440 


Transmitter 4 (Tr4) 


Saber 



Note that Tr2 and Tr4 are different radios, but they are of the same model. 
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(d) 



Figure 10: Preprocessed signal samples from four different transmitters: 
(a) Signal 1 from Transmitter 1, (b) Signal 1 from Transmitter 2, (c) Signal 
1 from Transmitter 3, (d) Signal 1 from Transmitter 4. 
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IV. PROCESSING RESULTS 



A. WAVELET TRANSFORM TECHNIQUE 

The data was re-processed by the wavelet transform technique to validate the 
integrity of the acquired data and the results of the previous work [8]. The distance 
measure algorithm is applied to classify or identify the signals (i.e., which signal 
belongs to which transmitter). The eighth-order Daubechies wavelet function is used 
to compute the wavelet transform of the processed signals [8]. The signal-to-noise 
ratios of all the signals are computed. The average SNRs are 34.8 dB, 41.1 dB, 40.6 
dB, and 34.6 dB for set 1, set 2, set 3, and set 4, respectively. Scale 11 appears 
to be the best scale for the data at hand [8]. Ten noisy realizations are generated 
using white Gaussian noise. The noise is added to the set to increase the noise level 
(i.e., to decrease the signal-to-noise ratio). The probability of identification versus 
different SNRs’ plots for the wavelet transform technique are obtained and shown in 
Figure 11 for the four transmitters. It can be seen that the WT distance algorithm 
is able to classify the four signal sets accurately. Also, good identification results are 
achieved at low signal to noise ratios. For example, it is found that for transmitter 1 
a 1.0 probability of identification is achieved at 16 dB, which is considered a reliable 
identification. For transmitters 2, 3, and 4, it is found that the WT distance algorithm 
is reliable above the SNR levels of 28 dB, 23 dB, and 21 dB, respectively. These results 
are identical to those obtained in [8]. 

B. NORMALIZED WAVELET TRANSFORM TECHNIQUE 

In this technique we scaled the signals and the templates of each transmitter to 
normalize the energy of the signals and the templates. The normalization is done by 
dividing the data vector ( mdmyi ) by the square root of its sum of the square values 
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(a) (b) 





(c) (d) 

Figure 11: Probability of identification (Pi) of the signals using the WT 
approach: (a) Pi of signals from Trl, (b) Pj of signals from Tr2, (c) Pj of 
signals from Tr3, (d) Pj of signals from Tr4. 
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as shown below 



mdmyii(n) 



mdmyii(n) 




,9 

,4 



(29) 



where mdmyit(n) is the i th pre-processed signal of the <? th transmitter. The chosen 
templates are normalized similarly. The wavelet transform is computed and the WT 
distance algorithm is applied. Ten realizations are conducted using different randomly 
generated noise for different SNRs. The results show that some of the WT distance 
between the template and its own set gets larger, while the WT distances to the other 
sets get smaller. This will not lead to a firm classification since a good and reliable 
classification can only be achieved if the distance between the template and its own set 
is small and its distance to the other sets are large. The results of using the normalized 
wavelet transform is shown in Figure 12. We see shows that this method performs 
better than the wavelet transform for two sets of signals and worse for the other two 
sets. For example, in the wavelet transform method a high probability of identifica- 
tion of transmitter 1 can be achieved at an SNR of 16 dB. While using the normalized 
wavelet transform method, a high probability of identification can be achieved at 21.8 
dB. Therefore, a degradation of 5.8 dB seems to occur. These results reveal that the 
normalized version of the wavelet transform does not perform better than that without 
normalization. This may be due to the negative values in some of transmitter 1 pre- 
processed signals. The negative values contribute to the square root of the sum of the 
squares. Signal 1 is selected as a template for the first set and, as shown in Figure 13, 
does not have large negative values. For illustration, Figure 13(b) shows signal 9 from 
set 1, where the negative spike can be observed. If we remove the negative spikes from 
the signals, the distance between the template and its own set become smaller and 
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(a) (b) 





(c) (d) 



Figure 12: Probability of identification (Pi) of the signals using normalized 
wavelet approach: (a) Pi of signal from Trl, (b) Pi of signal from Tr2, (c) 
Pi of signal from Tr3, (d) Pj of signal from Tr4. 



distances to other sets become larger, indicating that the negative spikes cause part 
of the performance degradation. 
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Figure 13: Plots of preprocessed signals 1 and 9 from Transmitter 1: (a) 
signal 1 from Trl, (b) signal 9 from Trl. 

C. SPECTRAL CROSS-CORRELATION TECHNIQUE 

The aim of this approach is to utilize an automated method to introduce a ro- 
bust identifier. The normalized cross correlation is bounded by +1 and —1. A “1” 
signifies that the two signals under consideration are 100% alike while “0” signifies 
no likeness. The normalized cross-correlation is independent of the signal ampli- 
tude. It measures the waveform similarities of two signals. We take the FFT of 
the preprocessed data, extract the magnitude squared of the coefficients, and retain 
the spectral coefficients with non-zero values. The dominant spectral coefficients are 
available from the first 32 spectral locations to capture the relevant information. The 
four averaged templates (Figure 7), chosen for their similarity within each set, are 
cross-correlated with the four signal sets. For the different SNRs, ten realizations are 
generated using independent noise. The processing results using this technique are 
shown in Figure 14 by plotting probability of identification (Pi) versus SNR. 

As we can see from the plots, the spectral correlation technique gives a better 
result than those based on the normalized and the non-normalized wavelet transform 
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Figure 14: Probability of identification (Pi) of the signals using spectral 
cross-correlation approach: (a) Pj of signals from Trl, (b) Pi of signals 
from Tr2, (c) P* of signals from Tr3, (d) Pi of signals from Tr4. 



(i.e., high probability of identification at lower signal-to-noise ratios) except for tem- 
plate 1 . The similarity between the waveform of template 1 and those of the signals 
belonging to transmitter 3 may be the reason behind the unsatisfactory performance 
of template 1. This similarity more likely occurs when we add white Gaussian noise 
to the signals of set 3. We investigated the data leading to the ROC curve of template 
1, which is shown in Figure 15. 
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Figure 15: Performance curve for the probability of proper classification 
(Ppc) versus the probability of miss-classification (P m c) of template 2 at 
19.8 dB. 

It was found that, to obtain a certain probability of classification, the thresh- 
old should be set to a specific value. For example, when the threshold is set to 
0.98 the probability of proper classification is found to be 0.3778 with a 0 proba- 
bility of misclassifi cation. When it is required to have a high probability of proper 
classification as 0.9778, the probability of mis classification is 0.0778 which can be 
considered as an acceptable case. This investigation was done with template 1 at a 



29 



signal-to-noise ratio of 19.8 dB. The other templates are performing perfectly, as 
shown in Figure 14, with low SNRs of 18.1 dB, 17.6 dB, and 11.6 dB for template 
2, 3, and 4, respectively. The threshold plots of the degraded signal-to-noise ratio of 
10 db and 22 dB for template 1 are shown in Appendix B (Figures 18 and 19). No 
further analysis is performed for templates 2, 3, and 4 because template 1 failed to 
work below 19.957 dB. 

D. PROCESSING SUMMARY 

The results of the three techniques are summarized in Figures 16 and 17. From 
the plots it can be seen that the spectral correlation approach outperforms the other 
two approaches for most transmitters. It was found that, for templates 2, 3, and 4, 
the spectral correlation gave a high probability of identification of 1.0 at low signal- 
to-noise ratio up to as low as 18.1 dB, 17.6 dB, and 11.6 dB, respectively. Whereas 
the wavelet and the normalized wavelet gave a high probability of identification of 1.0 
up to as low as 28.1 dB and 22.1 dB for template 2; 23.8 dB and 21.8 dB for template 
3; 20.6 dB and 24.6 dB for template 4. The spectral correlation approach works well 
for templates 2, 3, and 4. 
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(a) 




(b) 



Figure 16 : Probability of identification (Pj) plots versus signal-to-noise 
ratio (SNR in dB) for the three techniques, (a) Pi of signals from TRl, 
(b) P{ of signals from TR2. 



31 




(a) 




(b) 



Figure 17: Probability of identification (Pj) plots versus signal-to-noise 
ratio (SNR in dB) for the three techniques, (a) P\ of signals from TR3, 
(b) Pj of signals from TR4. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CLASSIFICATIONS OF PUSH-TO-TALK COMMUNICATIONS 
SIGNALS 

This thesis investigated the application of the wavelet transform and spectral 
correlation to identify several types of push-to-talk transmitters. The data provided 
for this work was turn-off /turn-on transients collected from push-to-talk radios. Nine 
transmissions from four different transmitters were used in the identification part of 
this thesis. The turn-on transients are unique for each transmitter. Results show that 
the spectral cross-correlation method out-performs the other two techniques. 

B. RECOMMENDATION FOR FUTURE STUDIES 

The spectral correlation approach introduced in Chapter IV, Section C gave 
the most promising results in classifying/identifying the four signal sets. An attempt 
was made to select a robust template from each signal set for wavelet processing by 
averaging the two dominant signals in each set with respect to the signal amplitude. 
However, poor results were observed. Further studies are recommended to find a 
robust template selection for wavelet processing to make the decision statistically 
more reliable. 

No attempt was made to use information from more than one scale simultane- 
ously. This might make the identification process more robust, so this subject is a 
goal for future study. Finally, other wavelet functions could be investigated for use 
in the identification process. 



33 




34 



APPENDIX A 



This program is to run wt signal to transform the data into a pulse shape form. It is 
also to map them and select the first four signals as templates. 



'/. DATARUN1.M 

'/, Must load data from thesis 1 subdirectory, 
'/, cd to ../thesis, then run datarunl 

'/, written by: abdulla 
'/, last mod: 02/28/96 

Ys=D; '/.define strings ’y -yyy’ 

Xs=D ; 



for n=l :4 

Xs=[Xs , ’x*] ; 

Ys= [Ys , ’y’]; 
dispC’ ’) 

dispC [’Outer loop: ’,Xs]) 
for m=l:9 

Sl=[Xs, int2str(m)] ; 

S2=[Ys, int2str(m)] ; 
disp([’ Inner loop: ’,S1]) 
evalC [* w* ,S1 , ’=mdm’ ,S2, ’ ; ; 
eval( [’a’ ,S1 , ’=map(w’ ,S1 , ’ ,8) ; ’] ) ; 
if m==l 

eval( [’ at ’ ,int2str(n) , ’=a’ , SI ; 

end 

end 

end 
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The purpose of this program is to compute the distance measurement between the signals 
and the four templates. 



*/, DATARUN2.M 

*/, THIS PROGRAM FIND THE DISTANCE BETWEEN THE SIGNALS AND THE TEMPLATES 

'/, written by Abdulla M. 

'/, last modified 02/28/96 



Xs= [] 

for n=l:4 

Xs=[Xs , ’x’] ; 
for m=l :9 

evalCC’dx* ,int2str(n) , * (* ,int2str(m) , * , : )=distance2(a’ ,Xs , int2str(m) , . . . 
’jat^,’)’]) '/, template variables are atl, at2, at3, at4. 

end 

end 
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This program transforms all noisy data into a pulse shape. 



'/. C0NV1.M 

'/, Must load datagen from thesisl subdirectory, 
'/, cd to ../thesis, then run datarunl 
'/, written by: abdulla, Abdulla M. 

'/, last mod: 05/4/96 

Ys=[]; '/.define strings ’y -yyy* 

Xs= □ ; 



for n=l :4 

Xs= [Xs , ’ x’] ; 

Ys= [Ys , ’y’]; 
disp( ’ ’) 

disp( [’Outer loop: ’,Xs]) 
for m=l:9 

Sl=[Xs, int2str(m)] ; 

S2=[Ys, int2str(m)] ; 
disp([’ Inner loop: ’,S1]) 
eval( [’mdm’ ,S2, ’=wtsig( ’ ,S1 ; 
end 

end 
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xc.m 

'/,*/, This program computes the fft of the pre-processed signals (mdmy’s) and 
cross correlate them with the chosen templates 
*/,'/, Written by Abdulla, Abdulla M. 

'/,*/, Last modified 06/7/96 

L=4096; '/, length of the vectors mdmy’s 2 to the power of N 

zll=tmpltl; '/, tmplt2, tmplt3, and tmplt4 
zll=zll-mean(zll) ; 

a2=mdmyl; */, mdmyyl, mdmyyyl, and mdmyyyyl. 

z2=abs (f f t (a2 ,L) . ~2) ; 

z22=z2( 1 :32) ; 

z22=z22-mean(z22) ; 

cl=xcorr(zll ,z22, ’coeff ’) ; 

a3=mdmy2 ; 

z3=abs(fft(a3,L) .~2) ; 
z33=z3(l : 32) ; 
z33=z33-mean(z33) ; 
c2=xcorr(zll ,z33, ’coeff’) ; 

a4=mdmy3 ; 

z4=abs(fft(a4,L) .~2) ; 
z44=z4(l:32) ; 
z44=z44-mean(z44) ; 
c3=xcorr(zll ,z44, ’coeff * ) ; 

a5=mdmy4 ; 

z5=abs (f ft (a5 ,L) . "2) ; 
z55=z5(l :32) ; 
z55=z55-mean(z55) ; 
c4=xcorr (zll ,z55 , ’ coeff * ) ; 

a6=mdmy5 ; 

z6=abs (f f t (a6 ,L) .“2) ; 
z66=z6( 1 : 32) ; 
z66=z66-mean(z66) ; 
c5=xcorr(zll ,z66, ’coeff ’ ) ; 

a7=mdmy6 ; 

z7=abs(fft(a7,L) . ~2) ; 
z77=z7 ( 1 : 32) ; 
z77=z77-mean(z77) ; 
c6=xcorr(zll ,z77 , ’coeff’) ; 
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a8=mdmy7 ; 

z8=abs (f ft (a8 ,L) . ~2) ; 
z88=z8 ( 1 : 32) ; 
z88=z88-mean(z88) ; 
c7=xcorr(zll ,z88 , ' coeff * ) ; 

a9-mdmy8 ; 

z9=abs(fft (a9 ,L) . ~2) ; 
z99=z9(l :32) ; 
z99=z99-mean(z99) ; 
c8=xcorr(zll ,z99 , ' coeff ’ ) ; 

al0=mdmy9 ; 

zlO=abs(fft(alO,L) .~2) ; 
zl010=zl0(l:32) ; 
zl010=zl010-mean(zl010) ; 
c9=xcorr (zll ,zl010 , * coeff * ) ; 
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APPENDIX B 



In this appendix two plots are shown. In Figure 18 we can see that at 10 dB degra- 
dation one can set a threshold between template 1 and the other transmitters (i.e., 
all other signals other than transmitter 1 are not similar to template 1). In Figure 19 
we can see that transmitter 3 signals are becoming similar to template 1 and starting 
to cross each other at 22 dB degradation. 




Figure 18: 10 dB degradation experimental automated threshold for tem- 
plate 1/set 1, set 2, set 3, and set 4. 
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Figure 19: 22 dB degradation experimental automated threshold for tem- 
plate 1/set 1, set 2, set 3, and set 4. 
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